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BACKGROUND 


Mississippi  State  University  (MSU),  in  collaboration  with  members  of  the  National  Science 
Foundation  Engineering  Research  Center  for  Computational  Field  Simulation  at  Mississippi  State 
University  (MSU)  has  established  a  collaborative  research  program  fimded  by  AFOSR,  McDonnell 
Aircraft  Company  (McAir),  and  McDonnell  Douglas  Research  &  Development  (MDRL)  and  Tele¬ 
dyne  Brown  Engineering  (TBE)  Corporation. 

The  broad-T>ased  objective  of  this  research  project  is  to  improve  understanding  of  fluid  phenomena 
associated  with  complex  three-dimensional  full  (e.  g.  complete  aircraft)  configurations  under  vari¬ 
ous  physical  environments.  In  particular,  the  research  emphasis  is  place  on: 

■  Grid  generation  associated  with  grids  about  complete  aircraft  configura¬ 
tions 

a.  Efficiency  and  Economy 

b.  Block/Zone  Strategies  and  Methodologies 

c.  Adaptive  Gridding 

■  Solution  algorithms 

a.  Block/Zone  Strategies 

b.  Solution  Adaptation 

c.  Parallel  Processing 

■  Training  of  C  O  applications  specialists 
a,  University/Industry  personnel 

The  key  element  of  this  project  is  the  technology  transfer  providing  one-on-one  interactions  be¬ 
tween  university  personnel  and  scientists  from  Vtirious  aerospace  industries.  This  program  is  a 
three-year  collaborative  effort.  This  report  summarizes  the  accomplishments  and  algorithms  de¬ 
veloped  during  the  execution  of  this  overall  broad-based  research  program. 

INTRODUCTION 

There  is  little  doubt  that  the  use  of  various  types  of  mathematical  modeling  has  significantly  en¬ 
hanced  almost  every  aspect  of  hardware  design,  including  not  only  evaluation  of  basic  concepts, 
but  also  in  the  prototyping  post-produchon  stages.  The  manner  in  which  modeling  affects  design 
processes  has  been  shown  to  reduce  cost,  risk,  and  enhance  the  overall  final  product.  A  particular 
type  of  modeling  which  has  increasingly  come  into  play  over  the  past  several  years  is  now  com¬ 
monly  referred  to  as  computational  fluid  dynamics  (CFD),  which  involves  numerically  solving  a 
set  of  partial  differential  equations  over  some  spatial  and  temporal  domailr  The  use  of  CFD  has 
been  shown  to  be  particularly  well-suited  to  the  aerospace  industry  (Ref.  1).  CFD  has  also  been 
extensively  utilized  in  the  ground  testing  of  many  aerospace  configurations  ranging  from  subscale 
aircraft  models  to  full-scale  missiles,  to  full-scale  turbine  and  rocket  engines  (Ref.  2-3). 

To  use  CFD  in  the  analysis  of  a  particular  fluid  flow  problem,  several  things  must  either  exist  or  be 
created  and  must  function  as  a  complete  system  to  be  useful:  (1).  generation  of  a  "computational 


grid",  where  the  spatial  region  of  interest  is  split  up,  or  discrehzed,  into  many  smaller  regions,  (2). 
execuhon  of  a  software  package  (sometimes  referred  to  as  a  "flow  solver"),  which  is  a  set  of  instruc¬ 
tions  making  up  the  numerical  algorithm  with  which  the  particular  set  of  partial  differential  equa¬ 
tions  are  solved,  (3)  execution  of  post-processing  software  designed  to  allow  the  rendering  of  large 
amounts  of  digital  information  into  graphical  or  analog  depictions  of  a  particular  flowfield  phe¬ 
nomena,  (4)  a  computer  powerful  enough  (in  terms  of  speed  and  memory  size)  to  execute  the  flow- 
field  solver  in  a  practical  amount  of  time,  and  (5)  persormel  skilled  in  the  use  of  this  software  and 
hardware  and  in  the  interpretation  of  computed  results. 

In  the  spirit  of  AFOSR's  initiative  to  promote  university /industry  collaborative  research,  the  proj¬ 
ect  entitled,  "Numerical  Flow  Simulation  for  Complete  Vehicle  Configurarion"  has  been  launched 
involving  McDonnell  Aircraft  Company,  McDonnell  Douglas  Research  and  Development  Compa¬ 
ny,  and  Teledyne  Brown  Engineering  Company  and  the  NSF  Engineering  Research  Center  fu* 
Computaaonal  Field  Simulation  at  Mississippi  State  University.  The  research  and  development  ac¬ 
tivities  of  this  effort  is  focused  on  the  advancement  of  full  Navier-Stokes  flow  simulation  algo¬ 
rithms  to  predict  more  accurately  and  efficiently  the  flow  phenomena  involving  complex  configu¬ 
rations.  In  particular,  the  research  emphasis  is  placed  on:  surface/ geometry  generation,  domain 
decomposition,  strategies,  solution  adaptive  grids,  grid  refinement,  parallel  processing,  and  nu¬ 
merical  algorithms  for  simulation  of  various  physical  environments  associated  with  complete  ve¬ 
hicle  configurations. 


OBJECTIVES 

The  broad-based  objective  of  this  research  project  is  to  provide  a  forum  through  which  university 
and  industrial  researchers  can  jointly  pursue  research  and  development  pertinent  to  the  simulation 
of  fluid  flowfields  (steady  and  unsteady)  about  complete  flight  vehicles  over  a  wide  range  of  flight 
conditions.  Upon  completion,  this  research  will  significantly  reduce  the  "response  time"  for  ad¬ 
dressing  complex  CAD  applications  and  will  increase  the  quality  of  simulahon  through  the  im¬ 
provement  of  numerical  algorithms,  surface  rendering  and  grid  generahon  schemes,  and  the  man¬ 
ner  in  which  graphical  techniques  are  used  to  interpret  computed  solutions.  Specifically,  research 
focused  upon  advancing  the  current  state  of  knowledge  regarding  the  following  items  is  being  con¬ 
ducted: 

a.  Grid  Generation 
Surface/ Geometry  Generation 
Grid  Optimization 

Domain  Decomposition  Strategies 
Adaptive  Gridding 

b.  Flowfield  Solution  Algorithms 
Parallel  Processing 

Block  Interface  Treatment 
Solution  Adaptive  Techniques 

c.  Training  of  CFD  Applications  Specialists 

d.  Applications;  The  present  approach  is  centered  around  performing 
computations  about  the  complete  F15  Figher  Aircraft. 


TECHNICAL  APPROACH 


The  primary  objective  is  to  improve  all  facets  of  the  CFD  technology  associated  with  efficient  and 
"optimal"  simulation  of  complete  flight  vehicle  configurations.  Numerical  geometry  grid  genera¬ 
tion  is  considered  the  most  labor  intensive  part  of  any  CFD  application.  In  fact,  at  present  it  cein 
take  significantly  more  labor  time  to  construct  a  CFD  grid  than  it  does  to  execute  the  flow  field  sim¬ 
ulation  code  on  the  grid  or  to  analyze  results.  Also,  80%-90%  of  the  grid  generation  labor  time  is 
usually  spent  on  geometry  processing.  Therefore,  initial  emphasis  is  placed  on  grid  generation  en- 
hancments  with  the  objective  of  increasing  efficiency,  quality,  and  produchvity.  To  this  end,  a  can¬ 
didate  complete  aircraft  configuration  was  selected.  Three  distinct  geometric  descriptions  of  the 
same  fighter  aircraft  with  distinct  refrence  and  scaling  criteria  were  obtained  in  the  form  of:  (i)  a 
CAD/CAM  geometry  description,  (ii)  a  set  of  forty-nine  sculptured  surfaces  and  (iii)  a  volume 
grid  around  a  wing/body/vertical-horizonal  tail.  These  configurations  were  studied  and  ana¬ 
lyzed  to  investigate  algorithm  development  which  could  expedite  the  grid  generation  process.  In 
this  respect,  graphical  user  interactions,  CAD  interface,  and  automatic/interacrive  multi-block  in¬ 
terface  construction  have  already  been  explored  in  the  general  purpose  grid  systems:  GRAPE  Vine, 
EAGLEView,  GRIDGEN,  ICEM,  IGG,  and  others.  However,  there  exists  a  need  to  develop  method¬ 
ologies  to  redistribute,  refine  and  automatically  remap  surface  grids  while  maintaining  the  fidelity 
of  the  geometrical  description  of  the  original  surfaces  along  with  robust  control  functions  for  ellip¬ 
tic  generation  and  grid  adaptive  schemes.  In  order  to  explore  all  research  issues  associated  with 
the  simulahons;  the  following  technology  development  must  be  undertaken: 

TA.  1 .  Surface  Grid  Generation 

Develop /adapt  methodologies  to  manipulate  and  define  surface  grids 
associated  with  the  complete  geometry.  These  surface  grids  must  satisfy  the 
requirements  associated  with  the  structured  physical  to  computational 
mapping  pertinent  to  domain  decomposition  strategies  under  consider¬ 
ation. 

TA.2.  Grid  Construction 

Develop/adapt  techniques  to  optimize  2D/3D  grids  associated  with  the 
complex  geometries  under  consideration.  This  optimization  must  result  in 
efficiency  and  economy  of  the  overall  grid  generation  process. 

TA.3.  Grid  Adaptation 

Develop/adapt  grid  adaptive  procedure  and  validate  with  respect  to  com¬ 
plex  configurations. 

TA.4.  Flowtield  Simulations 

Adapt  flowfield  solution  codes  and  simulate  flowfields  involving  various 
physical  environments.  This  must  include  increasingly  complex  configura¬ 
tions  such  as  F15e  and  cruise  missile  geometries. 

TA.5.  Domain  Decomposition  ^ 

Develop  domain  decomposition  strategies.  This  is  accomplished  as  fol¬ 
lows: 

a.  Generate  a  multiblock  grid  around  a  configuration  under  consider¬ 
ation.  Identify  block  boundaries/ surfaces  strictly  on  the  basis  of  geo¬ 
metrical  characteristics. 
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b.  Optimize  boundaries/ surfaces  mesh  for  improved  orthogonality  eind 
smoothness  and  hence  refine  the  grid  using  the  grid  optimization  pro¬ 
cess. 

c.  Simulate  multib’.t  ck  flowfield  utilizing  the  blocked  grid. 

d.  Analyze  and  compare  these  solutions  in  order  to  evaluate  "optimal" 
placemeiii  of  block  boundaries/surfaces.  This  willbe  accompli^ed  on 
the  basis  of  the  equidistribution  law. 

e.  P  ^.-generate  optimized  grid  and  repeat  steps  e  through  d  until  optimal  place¬ 
ment  of  boundaries/surfaces  is  achieved. 

T/..6.  Parallel  Processing 

Develop  parallel  processing  methodologies  for  the  single  block  and  multi¬ 
block  environment  being  studied. 

TA.7.  Validation  . 

Validate  solution  results  in  comparison  with  available  test  data  and  analyze 
the  pay-off  from  the  enhanced  technology. 

The  computer  programs  EAGLEView  (Ref.  4-6)  and  GENIE  (Ref.  7-9)  are  utilized  for  accomplish¬ 
ing  grid  generation  and  the  codes  PARC3D  (Ref.  10)  and  EAGLE  (Ref.  11-12)  are  used  for  flowfield 
simulations.  PARC3D  is  based  on  the  finite  difference  approxinruite  factorization  scheme,  however, 
EAGLE  is  based  on  finite  volume  Roe's  upwind  scheme. 

All  research  activities  should  be  performed  collaboratively  between  MSU  researchers  and  indus¬ 
trial  partners. 


ACCOMPLISHMENTS 

First  Year’s  Accomplishments 

During  the  first  year,  the  F15e  geometry  was  acquired  from  Wright  Patterson  AFB  and  General  Dy¬ 
namics.  The  grid  generation  and  flowfield  simulations  on  a  single  block  configuration  involving 
wing /fuselage  with  inlet  merged  with  fuselage  (no  flow  through  inlet)  were  performed.  These 
simulations  were  performed  with  M»  =  1.5,  a  =  0®  utilizing  finite  difference  and  finite  volume 
codes.  The  following  enhancements  were  made  in  grid  technology: 

Surfoce  Generation 

The  sculptured  cubic  B-spline  surface  representation  in  direct  and  inverse 
forms  Wcis  developed  for  surface  grid  generation. 

Grid  Optimization 

A  grid  optimization  based  on  the  best  characteristics  of  algebraic  and  ellip¬ 
tic  generation  system  was  developed  and  coded  for  the  two-dimensional 
applications.  f- 

Grid  Adaptation 

A  two-dimer\sional  grid  adaptative  algorithm  based  on  the  algebraic  tech¬ 
nique  and  quick  elliptic  optimization  algorithm  was  initiated. 

Zone  Subset  Grid 


II 
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The  EAGLEView  grid  generation  system  was  enhanced  for  subset  elliptic 
grid  generation.  This  was  accomplished  by  allowing  higher  order  continu¬ 
ity  at  the  subset  boundary  interface. 


Second  Year’s  Accomplishments 

During  the  second  year,  the  following  tasks  were  performed: 

Surface  Generation 

The  sculptured  cubic  B-spline  surface  representation  in  direct  and  inverse 
form  is  extended  to  the  general  Non-Uniform  Rational  B-spline  representa¬ 
tion. 

Grid  Optimization 

A  Hermite  transfmite  interpolation  algorithm  with  automatic  evaluatioirs 
of  slopes  and  proper  twist  vectors  is  developed.  A  grid  optimization 
algorithm  developed  during  first  year  is  extended  to  three-dimensional 
applications.  A  novel  approach  to  evaluate  control  functions  for  elliptic 
generation  system  is  developed. 

Grid  Adaptation 

A  two-dimensional  grid  adaptive  algorithm  is  developed  based  on  the 
equidistribuhon  principle  utilizing  the  best  characteristics  of  algebraic  or 
elliptic  general  system.  A  weight  function  is  formulated  using  the  Boolean 
sum  of  using  various  physical  characteristics  and  adaptation.  A  code  is 
completely  vadidated  and  trar^sferred  on  to  Teledyne  Brown  Engineering. 

GENIE  Enhancements 

In  order  to  interpret,  manipulate  and  define  the  F15e  CAD  geometry,  vari¬ 
ous  curve/surface  manipulation  routines  are  developed  in  GENIE.  These 
routines  are  added  to  the  development  of  CAGI  (Ref.  13)  code. 

Domain  Decomposition 

To  evciluate  criteria  for  domain  decomposition  strategies  a  grid-flowfield 
simulation  around  an  isolated  wing  is  performed. 

Interpolation/Search  Algorithm 

An  interpolation  code  for  post-processing  comparison  and  analysis  is  de¬ 
veloped. 


Third  Year’s  Accomplishments 

Redistribution  and  Remapping  Algorithms 

The  NURBS  representation  for  surfaces  is  utilized  to  redistribute  geometrically  complex  surfaces 
with  desired  point  distribution  criteria.  An  algorithm  to  automatically  remap  surface  grids  is  de¬ 
veloped.  The  remapping  algorithm  allows  the  automatic  generation  of  surface  grids  containing 
trimmed  subsurfaces,  gaps,  and  independent  grids  (e.g.  wing-body  intersection  on  the  body).  The 
body  surface  grid  can  be  automatically  generated  by  blending  the  intersected  wing  part  as  the  part 
of  the  grid. 
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Grid  Optimization 

The  surface  grid  optimization  based  on  the  Hermite  interpolation  parametric  space  associated  with 
NURBS  surface  is  developed.  The  elliptic  generation  system  with  the  control  functions  developed 
during  the  second  year  have  been  extendi  for  the  surface  grid  refinement. 


Grid  Adaptation 

The  algebraic  grid  adaptive  scheme  has  been  extended  to  three  dimensions.  The  redistribution 
scheme  along  with  weight  function  based  on  a  Boolean  sum  of  flow  characteristics  has  been  en¬ 
hanced  for  surface-by-surface  grid  adaptation  in  the  desired  direction  (or  combination  of  direc¬ 
tions  e.g.  I,  J,  or  J,  K  or  I,  K,  or  I,  J,  K).  Various  computational  examples  of  practical  interest  were 
exercised  using  this  grid  adaptation  scheme. 


Complete  F15  Configurations  &  Domain  Decomposition  Strategies 

The  grid  generation  process  associated  with  the  complete  F15  aircraft  configuration  has  been  com¬ 
pleted.  The  work  on  the  simulation  cind  determination  of  "optimal"  placement  of  zonal  interfaces 
is  on-going  or  continuing. 


Applications 

Three  different  F15e  grid  configurations  have  been  utilized  to  establish  a  baseline  computational 
capability  and  to  evaluate  enhancements  made  in  the  grid /simulation  methodology.  These  config¬ 
urations  are: 

i.  Wing/fuselage:  Single  block  and  four-block  configuration  -  580,000 
grid  points.  (Figure  la-d). 

ii.  Vertical-horizontal  tail/wing/fuselage:  Five  block  configuration  (Fig¬ 
ure  2a-b)  using  c-grid  on  the  fuselage.  This  grid  was  acquired  from 
McDonnell  Aircraft  Company. 

iii.  Vertical-horizontal  tail/wing/fuselage:  Five  block  configuration  using 
H-O  grid  on  the  fuselage  (Figure  3). 

iv.  Complete  configuration  body /iniet/nozzle  (Figures  4-20). 

In  the  wing/ fuselage  configurations,  approximately  580,000  grid  points  were  used  for  single  block 
(173  X  42  X  80  grid).  Approximately  2.2  million  grid  points  were  used  in  tl^C-H  tail /wing /fuse¬ 
lage  configuration. 

The  purpose  of  this  research  was  to  produce  a  detailed  three-dimensional  volume  grid  of  the  full- 
configuration  F-15  EAGLE.  The  grid  was  constructured  using  the  grid  generation  package  GE- 
NIE'^+,  an  algebraic  grid  system  developed  by  Dr.  B.  K.  Soni  at  the  National  Science  Foundation 
Engineering  Research  Center  located  at  Mississippi  State  University. 
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The  original  geometry  definition  for  the  F-15  was  supplied  in  the  form  of  49  surfaces,  each  describ¬ 
ing  a  distinct  part  of  the  aircraft.  These  surfaces  varied  in  dimension  and  orientation  with  respect 
to  the  plane  of  symmetry.  Using  GENIE'*"*',  the  surfaces  were  refined  and  fit  together  in  order  to 
maintain  continuity  along  grid  lines.  The  refinement  was  accomplished  using  the  Akima  spline 
routine  found  in  GENIE'*"*'.  This  spline  scheme  allowed  for  manipulation  of  the  surfaces  while 
maintaining  the  integrity  of  the  geometry  definition.  Once  the  grid  of  the  surface  geometry  was 
completed  and  consider^  of  acceptable  quality,  development  of  the  volume  grid  was  begun. 

The  first  surface  grid  was  a  wing/body /tail  combination.  This  grid  was  a  highly  accurate  represen¬ 
tation  of  the  fuselage,  wing,  and  vertical  cind  horizontal  tails,  but  did  not  attempt  to  model  the  inlet. 
Instead,  the  fuselage  was  smoothly  molded  into  the  wing,  simplifying  the  geometry.  This  grid  con¬ 
tained  5  blocks  and  was  a  C-type  grid.  Numerous  solutions  were  performed  on  this  grid,  and  one 
of  these  is  included  in  the  accompanying  figures. 

The  size  of  the  final  volume  grid  was  81  x  141  60,  or  685, 260  grid  points,  with  the  i-direction  wrap¬ 
ping  around  the  body,  the  f-direction  streamwise,  and  the  k-direction  normal  to  the  surface.  The 
surface  grid  included  the  hiselage,  the  inlet,  the  nozzle,  and  a  flow-through  chamber  cormecting 
the  inlet  to  the  nozzle.  In  order  to  properly  execute  a  solution  algorithm  on  the  grid,  the  volume 
grid  was  broken  into  5  blocks.  Accurate  modeling  of  the  inlet  and  nozzle  regions  was  considered 
of  high  priority  for  this  grid,  and  considerable  effort  was  invested  to  achieve  this  goal.  Numerous 
pictures  of  the  final  surface  and  volume  grids  are  included. 

Future  plans  include  completing  the  present  surface  grid  by  adding  the  wing  and  the  horizontal 
and  vertical  tails.  The  domain  decomposition  of  this  grid  may  prove  to  be  a  difficult  process  due 
to  the  complexity  of  the  geometry.  Also,  accurate  simulation  of  the  flow  through  the  inlet/ nozzle 
combination  will  be  a  strong  test  for  the  solution  algorithm. 

A  comparative  study  is  made  between  a  finite  difference  versus  a  finite  volume  code.  The  afore¬ 
mentioned  configuration  was  simulated  on  the  parallel  processor  in  collaboration  with  BRL  scien¬ 
tists. 


University-Industry  Collaboration 

The  broad-based  objective  of  the  present  collaborative  project  is  to  provide  a  forum  through  which 
university  and  industrial  researchers  can  jointly  pursue  research  and  development  pertinent  to  the 
simulation  of  fluid  flowfields  about  complete  flight  vehicles  over  a  wide  range  of  flight  conditions. 
To  this  end,  during  this  year,  the  following  activities  were  performed: 

•  Paul  Adams,  a  MSU  graduate  student  in  Aerospace  Engineering  spend  three 
semesters  (one  semester  during  the  first  year  and  two  semesters  during  the 
second  year)  at  McDonnell  Douglas  Research  and  Development  Laboratory 
(MDRL)  as  an  intern.  He  graduated  in  the  Fall  of  1992  with  a  master's  thesis 
entitled.  Decomposition  Strategies  for  A  Muitiblock  Time  Dependent  Grid 
About  A  Helicopter  Rotor/Fuselage  Configuration.  His  work  wassupported 
and  funded  by  MDRL.  The  experience  gained  in  this  study  will  beapplied  to 
complete  F 1 5  simulations. 

■  Paul  Craft,  an  MSU  graduate  student  in  Aerospace  Engineering  spent  a  se¬ 
mester  at  MDRL  as  a  summer  intern.  He  is  working  on  the  Domain  Decomposi¬ 
tion  Strategies  and  Interpolation/Search  techniques  applicable  to  multi¬ 
block  configurations.  He  will  present  his  Master's  thesis  in  December  1993. 


•  Various  trips  to  St.  Louis  were  mode  by  MSU  researchers.  In  particular,  five- 
three  day  trips  (on  average)  per  year,  to  St.  Louis  were  made  by  MSU  re¬ 
searchers.  This  has  facilitated  the  close  working  relationship  between  MSU, 

McDonnell  Aircraft  Company  (McAir),  and  MDRL  engineers. 

Two  graduate  students  and  one  undergraduate  student  visited  St.  Louis  with 
faculty  members.  Mr.  Bruce  Vu,  who  Isa  Ph.D.  student,  is  pursuing  his  research 
in  Domain  Decomposition  Strategies  on  F15e  configurations. 

■  Three  trips,  per  year,  to  Huntsville,  Alabama,  were  made  by  MSU  faculty  and 
graduate  students.  Mr.  Jiann-Cherng  Yang  has  completed  a  Ph.D.  in  com¬ 
putational  engineering  with  research  emphasis  on  grid  adaptation.  MSU  re¬ 
searchers  worked  very  closely  with  researchers  at  Teledyne  Brown  Engineer¬ 
ing  regarding  their  requirements. 

■  MDRL  McAir  and  Teledyne  Brown  Engineering  researchers  made  two  trips  to 
MSU.  These  trips  were  organized  in  conjunction  with  NSF  Engineering  Re¬ 
search  Center  for  Computational  Field  Simulation  industrial  programs  and 
collaborative  project  related  schedule. 

•  Two  trips  to  the  U.  S.  Army  Ballistic  Research  Laboratory  (BRL)  were  made  by 
MSU  researchers.  A  coilaborative  effort  has  progressed  between  BRL  and 
MSU  in  the  area  of  parallel  processing.  Two  trips  were  made  by  BRL  research¬ 
ers  to  MSU. 

■  The  work  associated  with  grid  quality  analysis  was  conducted  due  to  con¬ 
tractual  arrangement  between  McAir  and  MSU.  This  analysis  has  been  uti¬ 
lized  in  grid  generation  associated  with  multiblock  strategies  for  FI 5  com¬ 
plete  configuration. 

•  The  work  on  helicopter  configurations  has  been  conducted  due  to  con¬ 
tracted  arrangements  between  MDRL  and  MSU.  The  challenge  of  gridding 
helicopter  blade  rotation  in  a  multiblock  environment  has  made  a  singificant 
contribution  in  the  development  of  grid  generation  techniques  associated 
with  rotating  machinery. 

■  The  work  on  adaptive  gridding  is  conducted  due  to  contractual  arrange¬ 
ment  between  Teledyne  Brown  Engineering  and  MSU.  The  adaptive  algo¬ 
rithm,  especially  in  three  dimensions,  will  be  directly  applicable  to  FI 5  flow 
field  simulation  and  domain  decomposition  strategies. 

The  university-industry  collaboration  is  working  out  very  well  for  both  industry  and  the  universi¬ 
ty.  The  facility  at  the  Engineering  Research  Center  at  Mississippi  State  University  provides  an  ex¬ 
cellent  environment  for  industrial  visitors  as  well  as  for  this  overall  collaboration.  This  collabora¬ 
tion  has  been  proven  to  be  very  beneficial  to  participating  students.  Three  Ph.D.  students,  two 
masters  students,  and  four  undergraduate  students  have  participated  in  this  overall  effort. 


THEORETICAL  DEVELOPMENT  ^ 

NURBS  Surface  Representation 

The  Non-Uniform  Rational  B-spline  curve  and  surface  representation  (Ref.  14-15)  is  developed  in 
direct  and  inverse  formulation  for  redistributing  complex  surface  grids.  In  general  a  NURBS  curve 
is  expressed  as 
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(1) 


Y^WH)  P(i)b,it) 

am  =  .=2 - 

X  (0 


where  the  notation  is  interpreted  as  follows: 


The  IV (i)  are  the  weights  (positive  real  numbers). 
The  P(i)  are  the  control  points  ( points  in  ) . 
The  are  the  B  —  spline  basis  functions. 


A  NURBS  surface  (3-D)  can  be  expressed  parametrically  in  the  form  ( Ref.  14): 


where 


G  (s,t) 


K\  K2 

1=0 7=0 
ATI  /t2 

Y  Y^W(i,j)bi(s)bj(t) 

1  =  0  7  =  0 


(2) 


The  W(i,j )  are  the  weights  (positive  real  numbers). 
The  P(i,  j  )  are  the  control  points  (points  in  ). 


The  bi  are  the  B-spline  basis  functions  of  degree  Ml  determined  by  the  know  sequence  S  (-Ml), . 
. .,  S  (N1  +  Ml).  The  bj  are  the  B-spline  basis  functions  of  degree  Ml  determined  by  the  know  se¬ 
quence  T  (-M2),  ...,T(N2  +  M2).  Here,  N1  =  K1  -Ml  +  2  and  N2  =  K2  -M2  +  1. 

The  surface  itself  is  parametrized  for  U(0)  <  s  <  U(l)  and  for  V(0)  <  f  <  V(l)  where  S(0)  <  U(0) 
<  U(l)  <  S(N1)  andT(O)  <  V(0)  <  V(l)  <  T(N2). 

For  a  given  curve  with  a  set  of  data  points:  Xq,  xi, . . . ,  xl,  we  can  find  the  control  points  as  follows 
in  (Ref.  14): 


Define  the  chord  length  parameterization  as 


f- 


Aj  ^  II  Zljc,  ll 
^,+  1  IM^i+i" 


where  A I  =  Jc,  +  i  “ 


The  control  points  do,  di, . . . ,  dL  can  be  obtained  by  solving  the  matrix 
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where 

/•o  = 

n  =  +  J,. )  X,, 

''l  ~  ^3L-l' 
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(4) 


g  ^  (-^1-2  +  ^i-1  )  +  ^.-1  +  ^/+1  )  ^ 

•  ‘  zl,_2  +  A^_^  +  J,-  +  J,-  + 

4l,_i  +  Ai  +  Zl,+  , 

t 

and  then  set 


t 

t 


d_  1  —  Xq  ,  di^_^  —  Xi 

The  control  points  of  the  surface  are  determined  using  the  tensor  product  formula  associated  with 
both  V,  s,  and  f  parametrics.  . 


Algebraic  Surface  Grid  Optimization 

Let  r  =  (x,  y,  z)  be  the  physical  coordinates  aind  (|,  r\)  denote  the  computational  space  for  the  surface. 
Let  r  (s,  t)  =  (x(s,  t),  y(s,  t),  z(s,  t))  be  the  NURBS  parametric  representation  of  the  surface  with  (s,  t) 
as  the  parametric  space.  Now  as  demonstration  in  Uie  Figure  21,  there  exists  a  one  to  one  relation- 
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ship  between  physical  space,  parametric  space  cind  computational  space.  Using  this  argument  a 
surface  grid  optimization  algorithm  is  developed  in  view  of  optimizing  orthogonality. 

For  orthogonality 


r|  •  ffi  =  0  is  required 


i.e. 


( Tj  5^  +  r,  )  •  {rs  Sr,  +  rj  trj)  =  0 


i.e. 


■  rs)  s^  Sr,  +  (rs  '  r,)  (s^  tr,  +  Sr,)  +  (r,  ■  r,)  tr,  =  0 


(6) 


and  also 


represents  the  desired  distance  space  in  r\  direction 


i.e. 


(rs  s,,  +  r,  tr,)  •  ( r,  5,,  +  r,  tr, ) 


(7) 


i.e. 


(rs-rs)s^  +  2(rs  •  rt)  (Sr,  t,,)  +  r,  •  r,  (8) 

Now  the  parametric  space  is  fomi.Sa  ed  as  the  normalized  arc  length  as  presented  in  Figure  22. 
On  the  boundary  t)  =  1  and  t]  =  rimax/ 1|  =  0  and  therefore  equation  (6)  results  in 

(rs  •  rs)  s^  Sr,  +  {rs  •  r,)  s^  tr,  =  0  (9) 

Since  S|  can  be  computed  using  finite  difference  formulations  on  the  bovmdaries  t)  =  1  and  ti  =  ilmax/ 
the  equations  (8)  and  (9)  can  be  solved  for  s,,  and  t,,.  ^ 

Similarily,  the  equations 

can  be  solved  for  S|  and  t|  on  ^  =  1  and  ^  =  ^max  r-c^.daries.  This  information  can  then  be  utilized 
to  complete  Hermite  transfinite  interpolation  on  the  parametric  space  and  the  surface  grid  points 
can  be  re-evaluated  on  the  redistributed  parametric  space. 


t 


t 


and 


2 

r|  •  r|  =  {desired  distance  t]  ^  direction) 


The  computational  example  of  this  technique  is  presented  in  the  Figure  23  a-f .  As  it  can  be  seen 
in  the  Figure  23  the  optimized  surface  grid  is  orthogonal  near  boundaries. 


Robust  Hermite  Transfinite  Interpolation  Scheme 


t 

t 

t 


The  Transfinite  Interpolation  Method  (Ref.  4-9, 16)  is  a  widely  used  grid  generation  technique.  It 
is  based  on  a  Boolean  sum  of  the  one  dimensional  parametric  interpolation  schemes.  The  utiliza¬ 
tion  of  Hermite  cubic  representation  in  grid  generation  is  attractive  because  of  the  control  over 
slopes  near  boundary /boundary  surface.  However,  the  evaluation  of  these  slopes  is  extremely  cru¬ 
cial.  The  quality  of  the  resulting  grid  is  highly  influenced  by  the  magnitude  of  the  slopes  and  the 
twist  vectors  (cross  derivatives).  An  algorithm  to  evaluate  slopes  and  twist  vectors  is  developed. 
Let  r  =  (x,y,z)  and  Q  =  (|,  q,  and  5)  denote  the  physical  and  computational  space  respectively.  Then 
for  Hermite  interpolation  in  ^  direction  the  vector  should  be  obtained  at  ^  =  ^min  arid  ^  =  ^niax 
surfaces.  The  following  equations  are  utilized  in  evaluation  of  r^-. 


=  0 


(10) 


•  {Lr,  X  r^)  =  V 


#  Where  the  desired  cell  volume  is  calculated  from  the  initial  linear  TFI  grid.  The  twist  vectors  are 

calculated  using  the  following  equations. 
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ELLIPTIC  REFINEMENT 

Let  r  =  (xi,X2,  X3)  and  Q  =  ,  ^2/  denote  physical  and  computational  space  respectively.  Define 

covariant  a  *  1  contravariant  vectors  (Ref.  as  follows: 

a,-  =  covariant  vectors  (*  =  li2,3); 

£*  =  contravariant  vectors  =  1,2,3) 

gij  =  Q.i  *  %  =  gji  (i  =  1,2,3),  O'  =  1,2,3) 
gU  =  a'-  =  gi‘  0  =  1,2,3),  0  =  1.2,3) 

,  ^  =  detl  gij  I  =  [fl,  •  (£2  ^3)]^ 


(^y)|*  ~  (derivative  of  g^j  with  respect  to  ^  * ) 


—  ■  L^j  +  L^i  *  L^j^k 


i  =  1,2,3;  j  =  1,2,3;  and*  =  1,2,3. 


Using  (13),  the  following  statement  can  be  obtained 


_  ^  St  gij  )^*  "1”  ^gj 


i  =  1,2,3;  j  =  1,2,3;  k  =  1,2,3; 
Consider  a  three  dimensional  elliptic  grid  generation  system 


3  3 


E  +  E^‘i{‘  =  ® 


where 


«  =  l  7=1 


g  g  ^  gj  m  gkm  gj  n  gkn^ 

i  =  1,2,3;  j  =  1,2,3 


(/,  y,  k)  and  (f ,  m,  n)  cyclic 
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Our  aiialysis  for  evaluation  of  is  as  follows: 

In  order  to  evaluate  the  forcing  functions  k  —  1,2, 3,  usual  practice  is  to  write  (15)  as 


3 


*=1 


(16) 


q  =  1.2.3 

Using  (13),  Equation  (15)  can  be  rewritten  as  follows: 

3  3  3 

1  =  1  j=l  *=1 

3  3  (17) 

"  Z  Z  ® 

1=1  j=l 

q  =  1.2,3 

Observe  that,  using  (13)-<15),  Equation  (17)  can  be  written  in  terms  of  metric  terms  and  their  deriva¬ 
tives.  Also 


8ii  =  •  £^,  =  II  r^iV 

represent  an  increment  of  arc  length  on  a  coordinate  line  along  lines  which  varies  and 

8ij  =  £|.  •  £5>=||/-|i|l  •  ll'■^>ll  •  cose,i^j 

represent  measure  of  orthogonality  between  grid  lines  along  which  and  varies.  These  quanti¬ 
ties  can  be  evaluated  if  the  desired  increment  in  the  arc  lengdi  and  the  desired  angles  between  grid 
lines  are  known. 

Looking  at  the  "precise  control  of  spacing"  property  of  the  algebraic  grid  (Ref.  49),  the  quantities 
"gij"  can  be  evaluated  from  the  weU-defined  algebraic  grid,  and  using 

8ij  =  (  )  (y^  )  (cos  d)  ^ 

where  6  is  the  desired  angle  between  ^  grid  lines,  the  quantihes  g^j,  i  ^  j  can  be  evaluated.  Once 
all  g,y's  are  known,  the  Equation  (17)  can  be  solved  for  the  forcing  functions  0^,  k=  1,2,3.  In  par¬ 
ticular  when  orthogonality  is  assumed  i.e.,  6  =  90°  or  gy  =  0  for  i  j,  then  0^,  k  =  1, 2, 3  can  be 
forrr.ulated  as 
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1  -d_ 

2 


(in  j'  cyclic 

\  oii  SjjJ 


(18) 


A 


k=  1,2,3 

where 

ik  =  -  )  c>’c/ic 

dii  dyj 

i.j.k=  1,2,3 


(19) 


we  assume  that  gj^  ^  0  for  i  =  1, 2, 3 

An  application  of  these  cot>trol  functions  result  in  a  smooth/nearly  orthogonal  grid  in  fewer  itera¬ 
tions  (3-5)  of  the  elliptic  solver.  Two  examples  demonstrating  the  application  of  these  control  hmc- 
tions  (3  iterations)  on  surfaces  involving  smooth  and  sharp  concave  and  convex  regions  are  pres¬ 
ented  in  Figure  24a-f  and  Figure  25a-d.  It  can  be  observed  that  the  elliptic  grid  is  smooth  and  near 
orthogonal  in  both  cases. 

These  control  fimctions  <ire  applied  in  surface/ volume  grid  refinement.  The  surface  elliptic  system 
is  obtained  by  transforming  Eq.  (15)  into  the  parametric  form  (Ref.  49).  The  twD-dimensioi\al  para¬ 
metric  space  (distribution  mesh)  is  refined  and  dien  surface  coordinates  are  evaluated  from  NURBS 
representation.  The  detailed  mathematical  description  for  surface  grid  optimization  is  provided 
as  follows: 

QUASI-TWO-DIMENSIONAL  ELLIPTIC  SYSTEM 


Assumptions 

1 .  The  surface  is  represented  as  a  coordinate  surface 

C  =  constant 


2.  The  ^  and  lines  are  perpendicular  to  the  g  lines 

L  e.  g,3  =  0  and  g23  =  0 


3.  The  principal  curvature  of  the  g  lines  vanish  on  the  surface 


S22  -  2gi2  +  g,,  {rr,r,  -  <p  Cr,)  =  X  \  X  r,j  \  (r^  X  rr,) 

where  ^ 

^22^1  g  •  -  2gj2r^^  '  {r^xr^)  +  g,,r;,^  •  (r^xrr,) 


By  the  chain  rule  for  differentiation. 
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where 


^2  ~ 

A  =  a'  Kss  -  2/3'  rst  +  y'  r,, 
a'  =  x,^  +  +  ^2 

/S'  =  XsXt  +  ys  yt  +  z,  z, 


y'  =  xP-  +  >>,2  +  zs^ 


\  rsX  r,  P 


{xs  yt  -  ys  xtf  +  (xs  z,  -  zs  x,)^  +  {y,zt  -  Zs  y,f 


We  generate  the  grids  on  a  curved  surface  by  the  following  four  steps: 


1 .  Generate  initial  grid  by  algebraic  method  and  obtain  the  corresponding  distribution  mesh  in 

the  parametric  (s,  t)  space.  ^ 

2.  Compute  the  control  functions  ^  and  xp. 

3.  Solve  the  quasi-two-dimensional  elliptic  system  to  obtain  the  new  distribution  mesh  in  the 
parametric  (s,  t)  space. 

4.  The  grids  on  the  parametric  (s,  t)  space  are  transformed  to  the  physical  (x,  y,  z)  space  using 
the  Non-Uniform  Rational  B^pline  (NURBS)  representation  of  the  surface. 
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Surface  Redistribution  &  Remapping 


The  NURBS  representation  is  used  for  the  standard  surface  description.  The  convex  hull,  local  sup¬ 
port,  and  variation  diminishing  properties  (Ref.  14-15)  of  B-spline  fvmctions  contribute  to  the  gen¬ 
eration  of  the  well-distributed  smooth  grid.  The  parametric  space  associated  with  NURBS  is  trans¬ 
formed  as  the  distribution  (Norm<ilized  arc  length  based)  mesh  (Ref.  31,  32).  The  geometrical 
entities  which  are  expressed  into  Non-NURBS  form  can  be  exactly  or  approximately  (with  desired 
tolerance)  converted  to  NURBS  presentarion.  The  evaluation  of  derivatives  r^.  r^.  r„.  Vff.  rff . . .  etc 
can  be  readily  from  the  formulation  (1).  The  redistributed  surface  grid  is  accomplished  by  evaluat¬ 
ing  the  NURBS  surface  at  the  respective  parameter  associated  with  the  desired  distribution  space. 
For  example,  the  network  of  control  points  and  the  associated  parametric  space  are  presented  in 
Figure  26a-b.  A  candidate  desired  distribuHon  mesh  is  demonstrated  in  Figure  27a.  The  resulting 
redistributed  surface  grid  is  presented  in  Figure  27b. 

It  is  important  to  note  that  the  redistributed  surface  grid  is  obtained  by  evaluating  the  NURBS  sur¬ 
face  at  the  desired  distribution  points.  Also,  there  exists  a  one-to-one  correspondence  between  the 
parametric  space  and  the  control  net.  The  straightforward  search  algorithm  is  utilized  to  evaluate 
a  cell  in  which  a  desired  parametric  point  resides.  A  redistributed  surface  resulting  from  a  regular 
cubic  spline  surface  is  presented  in  Figure  27c.  The  NURBS  surface  is  smooth  and  has  kept  the  fidel¬ 
ity  of  the  original  surface. 

Now,  consider  a  sculptured  surface  presented  in  Figure  28a.  Tne  surface  grid  is  mapped  as  a  C-grid 
on  the  surface.  However,  if  the  user  is  interested  in  mapping  an  O-type  grid  on  the  surface  then 
a  remapping  process  is  accomplished  by  essentially  creating  a  distribution  space  which  when  eva¬ 
luated  results  in  a  surface  grid  of  Figure  28b.  The  associated  parametric  (distribution)  spaces  are 
presented  in  Figures  29a-b.  The  creation  of  the  distribution  space  for  remapping  is  intuitive  and 
highly  application  dependent.  The  second  example  describing  the  process  of  remapping  an  O- 
type  grid  into  ah  H-t^e  grid  is  presented  in  Figure  30a-d.  Remapping  of  surface  grid  is  also  re¬ 
quired  when  the  interior  object(s)  is  (are)  to  be  kept  fixed  as  part  of  the  interior  surface  grid.  For 
example,  consider  the  control  net  and  the  assocaited  parametric  space  presented  in  Figure  a-b31. 
An  interior  object  on  the  surface  is  presented  in  Figure  32.  The  remapping  process  is  needed  to 
blend  an  interior  object  as  a  part  of  the  surface  grid.  An  interpolation  -  search  algorithm  based  on 
the  NURBS  presentation  is  developed  to  evaluate  parameters  associated  with  the  interior  object. 
This  evaluation  is  demonstrated  in  Figure  33.  The  interpolation/ search  algorithm  utilizes  deriva¬ 
tives  and  and  Taylor's  expansion  to  inversely  evaluate  parameters.  An  automatic 

algorithm  based  on  the  weighted  transfinite  interpolation  (Refs.  7-9)  in  two  dimensions  is  devel¬ 
oped  which  blends  the  parameters  associated  with  interior  object  into  an  overall  distribution  space. 
The  resulting  re-parameterized  distribution  space  is  demonstrated  in  Figure  34.  The  surface  grid 
is  then  evaluated  with  respect  to  this  new  distribuhon  space.  The  fidelity  of  the  geometry 
associated  with  the  interior  object  is  kept  precisely  on  the  surface  grid.  The  resulting  remapped 
surface  grid  is  presented  in  Figure  35. 


Algebraic  Surface  Grid  Optimization 

Detailed  mathematical  description  of  this  algorithm  is  provided  in  Appe?idix  A.4.  A  graphical 
view  of  linear  TFI  v/s  hermite  TFl  is  presented  in  Figure  36a-d  on  an  arbitrary  configuration.  The 
grid  lines  are  smoother  and  near  orthogonal  in  the  case  of  Hermite  TFl. 

These  formulations  are  coded  in  the  general  purpose  grid  system  GENIE.  Various  computational 
examples  have  been  exercised  to  demonstrate  the  success  of  this  technique.  These  examples  are 
provided  in  Appendix  A.8  in  a  formal  paper. 


Grid  Adaptation 


The  algorithm  and  a  computer  code  for  two  dimensional  grid  adaptation  is  completed.  This  algo¬ 
rithm  is  based  on  the  equidistribution  law  utilizing  NLTRBS-algebraic  technique  and  elliptic  gen¬ 
eration  system  with  the  forcing  fimctions  presented  in  previous  section.  This  work  is  supported 
by  Teledyne  Brown  Engineering  Corporation.  A  technical  AIAA  paper  containing  the  detailed 
mathematical  analysis  and  computational  examples  is  presented  in  Appendix  A.3  and  A.5.  The 
definition  of  weighting  function  using  the  linear  combination  of  Boolean  sum  of  various  flowfield 
characteristics,  NURBS-algebraic  regridding  process  maintaining  the  fidelity  of  the  associated  sur¬ 
face  geometry  and  elliptic  smoothing  and  rear  orthogonality  balance  are  the  salient  features  of  this 
algorithm. 

Grid  Adaptation 

With  structured  grids,  the  adaptive  strategy  based  on  redistribution  is  by  far  the  most  simple  to  im¬ 
plement,  requiring  only  the  regeneration  of  the  grid  and  interpolation  of  flow  properties  at  the  new 
grid  points  at  each  adaptive  stage  without  modification  of  the  flow  solver  unless  time  accuracy  is 
desired.  Time  accuracy  can  be  achieved,  as  far  as  the  grid  is  concerned  by  simply  transforming  the 
time  derivatives,  thus  adding  convective-like  terms  that  do  not  alter  the  basic  conservation  form 
of  the  PDEs. 

Adaptive  redistribution  of  points  traces  its  roots  to  the  principle  of  equidistributin  of  error  [50]  by 
which  a  point  distribution  is  set  so  as  to  make  the  product  of  the  spacing  and  a  weight  function 
constant  over  the  points: 


wJx  =  constant  (20) 

With  the  point  distribution  defined  by  a  function  |i,  where  ^  varies  by  a  unit  increment  between 
points,  the  equidistribution  principle  can  be  expressed  as 

wx^  —  constant  (21) 

This  one  dimensional  equation  can  be  applied  in  each  direction  in  an  alternating  fashion  [41].  A 
Direct  extension  to  multiple  dimensions  using  algebraic  [32],  variational,  and  elliptic  system  [33] 
has  been  developed. 

Weight  Function 

The  weight  function  is  a  very  important  part  in  the  overall  adaptive  process.  A  generalized  weight 
function  applicable  to  various  flow  field  characteristics  has  been  developed.  The  weights  are  com¬ 
puted  in  all  computational  directions  and  then  coupled  adaptation  is  applied.  A  linear  combina¬ 
tion: 


N 

1  -t-  ^j)  ^tf]  disf,  ^  i  ^ 

y=i 

where  N  =  number  of  flow  variables  (e.  g.  pressure,  temperature,  density,  etc.) 

kj  -  weighting  factor  associated  with  flow  parameter  >  0  (23) 


(Oj  =  ttj  qj  ©  kj  =  aj  qj  +  kj  +  (Uj  +  ^j-  D  kj 

qj  -  scaled  gradient  of  the  flow  variable  /  such  that  0  <  <  1 

kj  —  scaled  curvature  values  of  the  flow  variable  j  such  that  0  <  kj  <  I 
wtf-  weight  factor  that  enhances  the  total  effect  of  high  weighted  areas 
disf-  distribution  factor  that  can  keep  the  original  distribution 
0  <  <  1,  0  <  <  1. 


is  formulated  as  a  weight  function  utilizing  Boolean  sum  of  contributions  from  scaled  gradients 
emd  ctirvatures.  The  value  of  the  weight  contribution  is  controlled  by  the  weight  factors  and  is  at 
a  maximum  when  gradients  and/or  curvature  values  are  at  a  ntuiximum.  An  appropriate  scaling 
scheme  [32]  for  weight  factors  Oj,  Pi  has  been  developed  for  proper  weighting. 

Algebraic  Technique 

The  redistributed  algebraic  grid  is  generated  by  utilizing  a  surafce/ volume  distribution  mesh  as 
the  re-parametized  space  associated  with  Non-Uniform  Rational  B-spline  surface/volume  repre¬ 
sentation.  The  convex  hull,  local  support  and  variation  diminishing  properties  of  B-splines  make 
NURBS  an  excellent  choice  for  modeling  complex  geometries.  The  application  of  inverse  NURBS 
formulation  [6]  allows  reevaluation  of  control  points  which  influences  the  fidelity  of  solid  surface 
geometry  during  redistribution  process. 

Elliptic  Technique 

The  elliptic  generation  system: 


I  I 


1=1  j=\ 


+ 


=  0 


where 

r 

:  position  vector 

8'^ 

:  contravariant  metric  tensor 

:  curvilinear  coordinate 

Pk 

:  control  function 

(24) 


is  widely  utilized  for  grid  generation  [9] .  The  control  of  the  characteristics  and  distribution  of  a  grid 
system  Cem  be  achieved  by  varying  the  values  of  the  control  functions  Pjt  in  Equation  24  [6].  The 
application  of  the  one  dimensional  form  of  Equation  4  with  Equation  20  res|dts  in  the  definition  of 
control  functions  in  three  dimensions. 


These  control  functions  were  generalized  by  Eiseman  [55]  as 


(25a) 
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(25b) 
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In  order  to  conserve  the  geometrical  characteristics  of  the  existing  grid,  the  definition  of  control 
functions  is  extended  as 


(  ^initial  geometry  )  i  ^  )  (  *  1.2,3) 


(1)  (G)  (0) 

(  ^initial  geometry  )  ~  (  ^initial  geometry  )  *  (  ^t  ) 


(26) 


Grid  Characteristic  Analysis 

A  graphically  interactive  grid  characteristic  analysis  software  is  developed.  This  system  interac¬ 
tively  allows  surface  grid  to  be  manipulated  and  analyzed.  The  manipulation  operations  involve: 
picking  line,  reverse  ^  and  Tidirections,  remove  a  grid  line,  move  a  point,  and  perform  various  trans¬ 
formations.  The  quality  analysis  include:  aspect  ratio,  smoothing  measures,  orthogonality  mea¬ 
sures,  area  check  with  crossing,  and  measures  due  to  interaction  with  solution  features  and  grid. 
A  pictorial  view  of  angle  checking  is  provided  in  Figirre  37.  A  point  moving  operation  and  a  line 
removing  operation  is  demonstrated  in  Figures  38.  TTie  grid  manipulation  operation  displays  sur¬ 
face  geometry  (three-dimensional)  along  with  its  parametric  space  (two-dimensional)  as  demon¬ 
strated  in  Figures  39.  This  allows  extremely  easy  point  and  line  picking  due  to  the  two-dimensional 
parametric  space.  There  is  a  one-to-one  correspondence  between  a  point  on  the  three-dimensional 
surface  and  a  point  on  the  two-dimensional  parametric  space.  An  example  of  checking  surface 
crossing  is  presented  in  Figure  40.  The  list  of  crossings  found  is  provided  in  Figure  41.  This  list 
includes  (i,  j)  location  of  possible  cells. 

This  work  is  supported  by  a  McDonnell  Aircraft  Company. 


Domain  Decomposition  Strategies  &  Interpoiation/Search  Algorithms 

A  study  to  understand  placement  of  a  surface  interface  between  blocks  is  performed  considering 
flowfield  simulation  around  a  wing.  This  is  accomplished  by  grid  adaptive  algorithm  with 
NURBS  surface  re-distribution  technique  to  evaluate  the  placement  of  blocks  interfaces.  The  grid 
sizes  are  kept  small  in  order  to  simulate  flow  on  the  workstation.  Detailed  description  on  this  devel¬ 
opment  is  provided  in  Appendix  A.6.  A  program  to  interpret  flowfield  properties  at  a  cutting  plane 
is  developed.  The  program  uses  interpolation /search  techniques  to  evaluate  intersecting  points 
between  cutting  plane  and  the  volume  and  then  interpolate  flow  propertie^t  those  points.  A  de¬ 
tailed  description  of  this  development  can  be  found  in  Appendix  A.7.  Two  additional  technical  pa¬ 
pers  associated  with  grid  generation  are  provided  in  Appendix  A.8  and  A.9. 

FLOWFIELD  SIMULATIONS 

Flow  field  simulations  associated  with  the  aforestated  configurations  were  made.  These 
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i.  Wing/fuselage  -  single  block  configuration  (Figures  a-d) 

ii .  Vertical  -  horizontal  tail/wing/fuselage  -  five  block  configuration  (Fig¬ 
ures  2a-b)  using  c-grid  on  the  fuselage.  This  grid  was  acquired  from 
the  McDonnell  Aircraft  Company. 

iii.  Vertical  -  horizontal  tail/wing/frjselage  -  five  block  configuration  us¬ 
ing  H-O  grid  on  the  fuselage  (Figures  3). 

iv.  Complete  configuration  body /inlet/nozzle  (Figures  4-20). 

In  the  wing/fuselage  configurations,  approximately  580,000  (173  x  42  x  80  grid  for  single  block) 
grid  points  were  used.  Approximately  2.2  million  grid  points  were  used  in  the  C-H  tail /wing/  fu¬ 
selage  configuration. 

A  technique  to  re-distribute  surface  grids  using  Non-Uiuform  Rational  B-spline  (NURBS)  repre¬ 
sentation  has  been  developed.  The  NURBS  representation  allows  the  redistribution  of  grid  points 
on  the  surface  while  maintaining  fidelity  of  the  surface  geometry  (Ref.  31).  Flowfield  simulations 
using  similar  physical  conditions  and  grids  were  also  performed  on  the  Connection  Machine 
(CM-200)  in  parallel.  A  redistributed  surface  grid  for  the  vertical-horizontal  tail/ wing/ fuselage, 
configuration  iii,  is  presented  in  Figure  3. 

A  comparative  study  is  made  between  two  popular  algorithms. 

The  first  algorithm  is  based  on  Roe's  approximate  Riemann  solver  as  implemented  by  Whitfield 
(Ref.  20).  This  scheme  solves  the  time-dependent  Reynolds-averaged  Navier-Stokes  equations  in 
generalized  time-dependent  curvilinear  coordinate  systems.  Acciuacy  is  achieved  by  a  conserva¬ 
tive  finite-volume  discretization  which  satisfies  the  geometric  conservation  law  generalized  to 
moving  coordinate  systems.  The  flow  solver  code  is  written  in  a  block-structured  form  allowing 
calculations  with  nearly  unrestricted  arrangement  of  arbitrarily  sized  blocks  (Ref.  21).  Upwind  dif¬ 
ferencing  is  used  for  the  inviscid  terms  and  centered  differencing  is  used  for  the  viscous  and  diffu¬ 
sive  terms.  The  inviscid  flux  vectors  in  the  residual  are  computed  using  the  flux  difference  split 
theory.  To  summarize  the  approach.  Roe  averaging  (Ref.  22)  was  used  to  determine  the  numerical 
flux  at  cell  faces  with  the  schemes  of  Osher  and  Chakravarthy  (Ref.  23)  applied  to  achieve  second 
or  third-order  spatial  accuracy  .  There  is  some  deviation  in  construction  of  the  higher-order 
schemes  from  that  put  forth  in  (Ref.  23)  as  discussed  in  (Ref.  20).  In  keeping  with  the  upwind  nature 
of  evaluating  the  inviscid  flux  vectors  in  the  residual,  the  implicit  operator  is  also  upwinded.  How¬ 
ever,  it  was  found  in  (Ref.  20)  that  convergence  was  much  improved  if  flux  vector  splitting  (Ref.  24) 
was  used  in  the  implicit  operator  of  the  equation  rather  than  flux  difference  splitting.  The  viscous 
and  diffusive  terms  are  treated  explicitly,  and  the  turbulence  effects  are  modeled  with  the  well- 
known  Baldwin  and  Lomax  mixing  length  model  (Ref.  25),  although  present  compilations  were 
made  using  an  inviscid  model.  A  mathematical  description  of  this  algorithm  is  provided  in  Appen¬ 
dix  All. 

The  second  code  applied  is  the  PARC3D  code  as  developed  by  Cooper  and  Sirbaugh  (Ref.  26)  at 
the  Arnold  Engineering  Development  Center  (AEDC).  This  code  also  solves  either  the  Navier- 
Stokes  equations  or  the  Euler  equations,  subject  to  user  specification.  PARC3D  is  a  finite-difference 
code  which  applies  standard  central  difference  discretization  of  the  system  of  PDE's.  Jameson  style 
artificial  dissipation  is  applied  to  maintain  stability.  The  resulting  system  of  equations  was  solved 
herein  by  application  of  Pulliam's  (Ref.  27)  diagonalized  Beam  and  Warming  approximate  factor- 


ization  algorithm.  Among  other  characteristics,  PARC3D  utilizes  conservative  metric  differencing, 
offers  local  time-step  options,  and  provides  a  generalized  boundary  condition  treatment.  The  code 
can  be  applied  in  either  a  single  or  multi-block  mode.  This  results  in  a  code  which  is  quite  general 
and  easy  to  apply  to  complex  configurations. 

Since  the  architecture  of  the  Single  Instruction  Multiple  Data  Stream  Machine  (SIMD)  is  well  suited 
for  an  explicit  type  method,  an  explicit  multi-stage  multi-block  full  Navier-Stokes  solver  devel¬ 
oped  by  Patel,  Sturek,  and  Smith  (Ref.  28)  is  utilized  for  parallel  computations.  This  code  allows 
multi-zone  overlaid  grids.  To  maintain  stability  adaptive  artificial  dissipation  terms  based  on  a 
directional  eigenvalue  scaling  (Ref.  29)  have  been  added  to  the  equations.  A  local  time  stepping 
scheme  which  is  well  suited  for  pcirallel  processing  is  developed. 

Preliminary  computations  have  been  performed  simulating  the  flow  field  about  the  F15e  wing/fu¬ 
selage  configuration  at  a  Mach  number  of  1.5  and  an  angle  of  attack  of  0.0.  The  simulation  for  the 
grids  a  and  in  including  the  vertical  and  horizontal  tail  is  underway. 

The  simulations  were  made  by  application  of  the  previously  mentioned  Roe  finite  volume  scheme 
(EAGLE)  and  the  Beam-Warming  finite  difference  scheme  (PARC3D).  The  PARC3D  simulation 
was  made  as  a  single  block  computation  whereas  the  EAGLE  simulation  was  made  on  the  same 
grid  split  into  two  contiguous  blocks.  This  was  done  to  satisfy  in-core  memory  requirements  for 
the  UBIFLOW  code. 

Macroscopic  flow  features  for  the  two  computations  agree  quite  well  as  evidenced  by  comparing 
Mach  number  contours  on  the  fuselage  and  wing  surfaces  (Figures  42  and  43).  The  magnitude  of 
the  Mach  number  as  well  as  the  compression  and  expansion  patterns  match  reasonably  well.  Close 
examination  of  the  respective  contours  reveals  the  Roe  scheme  as  an  apparently  less  dissipative 
solution  which  is  consistent  with  anticipated  results.  Representative  quantitative  comparison  of 
computed  wing  pressure  distributions  are  shown  on  the  upper  surface  at  50%  chord  and  on  the 
lower  surface  at  25%  chord  in  Figures  44  and  45.  These,  and  other  supporting  figures,  illustrate  sim¬ 
ilar  overall  levels  of  agreement.  However,  disagreement  is  observed,  most  notably  on  the  lower 
surface  near  the  wing  tip.  Also  general  disagreement  is  noted  near  the  leading  and  trailing  edge. 

The  computations  at  a  Mach  number  of  0.98  and  an  angle  of  attack  of  1.1  degrees  will  be  compared 
with  experiment.  These  simulations  are  being  made  on  a  more  realistic  grid  (li  and  in)  representing 
the  fuselage/ wing/vertical  and  horizontal  tail  assembly.  Integrated  parameters,  such  as  lift  coeffi¬ 
cient,  will  also  be  compared  to  assess  the  overall  significance  of  discrepancies  within  the  solutions. 
These  computations  are  partially  complete  and  a  request  for  access  to  experimental  data  has  been 
made. 

An  additional  demonstration  of  capability  is  made  by  presenting  results  of  a  simulation  of  the  fuse¬ 
lage/wing/vertical  and  horizontal  tail  assembly  by  application  of  an  explicit  CFD  code  executing 
on  a  parallel  computer  architecture.  The  computation  was  made  to  simulate  conditions  at  Mach 
1.5  and  0.0  angle  of  attack.  Representative  results  are  shown  in  the  form  of  surface  pressure  dis¬ 
tribution  (Figure  46).  These  results  will  be  compared  with  computations  made  with  the  implicit 
Roe  scheme. 

One  results  of  this  research  was  the  development  of  various  advanced  tecl^iiques  needed  to  effi¬ 
ciently  solve  problems  encountered  during  the  F15  grid  generation  process. 

REFERENCES 

1.  Kutler,  R,  “A  Perspective  of  Theoretical  and  Applied  Computational  Fluid  Dynamics”, 
AIAA-83-0037,  1983. 


27 


2.  Phares,  W.  J.,  Cooper,  G.  K.,  Swafford,  T.  W.,  Jones,  R.  R.,  “Application  of  Computational 
Fluid  Dynamics  to  Test  Facility  and  Experimental  Design”,  AIAA-86-1733,  June  1986. 

3.  Droegemeier,  K.  K.,  “Numerical  Simulation  of  Thunderstorm  Outflows  and  Microbursts”, 
CRAY  CHANNELS,  p.  18,  Summer  1987. 


4.  Thompson,  J.  F.,  “Composite  Grid  Generation  for  General  Three-Dimensional  Regions  - 
The  EAGLE  Code",  AIAA  Journal,  Vol.  26.  No.  3,  p.  271,  March  1988. 

5.  Thompson,  J.  F.,  “A  Composite  Grid  Generation  for  General  Three-Dimensional  Regions  - 
The  EAGLE  Code",  AIAA  Journal.  Vol.  26.  No.  3.  p.  271.  March  1988. 


6.  Thompson.  J.  F.,  “Program  EAGLE-numerIcal  grid  Generation  System  User’s  Manual”, 
Vols.  II,  III,  AFATL-TR-87-15,  March  1987. 

7.  Soni,  B.  K.,  “GENIE;  GENeration  of  Computational  Geometry-Grids  for  Internal-External 
Flow  Configurations”.  Proceedings  of  the  Numerical  Grid  Generation  in  Computational  Fluid 
Dynamics^SB,  Miami, FL,  1988. 

8.  Soni,  B.  K.,  McClure,  M.  D.  and  Mastin,  C.  W.,  “Geometry  and  Generation  in  N+1  Easy 
Steps”,  The  First  International  Conference  on  Numerical  Grid  Generation  in  Computational 
Fluid  Dynamics  '68,  Miami,  FL,  1988. 

9.  Soni,  B.  K.,  “Two  and  Three-Dimensional  Grid  Generation  for  Internal  Flow  Applications  of 
Computational  Fluid  Dynamics”,  AIAA-85-1526. 

10.  Cooper,  G.  K.,  “The  PARC  Code:  Theory  and  Usage”,  AEDC-TR-87-24,  Arnold  AFB,  TN, 
October  1987. 


11.  Whitfield,  0.  L.,  “Implicit  Upwind  Finite  Volume  Scheme  for  the  Three-Dimensional  Euler 
Equations”,  Engineering  and  Industrial  Research  Station  Report,  Mississippi  State  Universi¬ 
ty,  MSSU-EIRS-ASE-85-1,  September  1985. 

12.  Whitfield,  D.  L.,  Janus,  J.  M.  and  Simpson,  L.  B.,  “Implicit  Finite  Volume  High  Resolution 
Wave-Split  Scheme  for  Solving  the  Unsteady  Three-Dimensional  Euler  and  Navier-Stokes 
Equations  on  Stationary  or  Dynamic  G  rids”.  Engineering  and  industrial  Research  Station  Re¬ 
port,  Mississippi  State  University,  No.  MSSU-^IRS-ASE-88-2,  August  1987. 

13.  Soni,  B.  K.,  Yu,  T.-Y,  and  Vaughn,  D.,  “CAGI:  Computer  Aided  Grid  Interface:  A  CFD  Pre- 
Processor",  Proceedings  of  the  CFD  Workshop  for  Rocket  Propulsion  Applications,  NASA 
Marshall  Space  Flight  Center,  Huntsville,  AL,  April  1992. 

14.  Farin,  G.,  “Curves  and  Surfaces  for  Computer  Aided  Geometric  Design:  A  Practical  Guide”, 
Academic  Press,  1990. 

15.  deBoor,  C.,  “A  Practical  Guide  to  Splines”,  Springer,  1978. 

16.  Steinbrenner,  J.  and  Foiits,  C.,  “GRIDGEN;  User’s  Manual”,  Wright  Research  and  Develop¬ 
ment  Center  and  General  Dynamics,  Fort  Worth  Division,  1 990. 

17.  Sorenson,  R.  L.,  “Three-Dimensional  Zonal  Grids  About  Arbitrary  Shapes  by  Poisson’s 

Equation”,  Numerical  Grid  Generation  for  Computational  Fluid  Mechanics  ‘88,  p.  88,  Piner- 
idge  Press.  y- 

1 8.  Steger,  J.  L.  and  Chaussee,  D.  S.,  “Generation  of  Body-Fitted  Coordinates  Using  Hyperbolic 
Partial  Differential  Equations",  SIAM  Journal  of  Scientific  Statistic  Computation,  p.  431 , 
1980. 

19.  Soni,  B.  K.,  Thompson,  J.  F.,  Stokes,  M.  and  Shih,  M.-H.,  “GENIE  ++,  EAGLEView  and  TI¬ 
GER:  A  General  Purpose  as  Special  Purpose  Graphically  Interactive  Grid  Systems", 
AIAA-92-0071. 


28 


20.  Whitfield,  D.  L.,  Janus,  J.  M.  and  Simpson,  L.  B.,  “Implicit  Finite  Volume  High  Resolution 
Wave-Split  Scheme  for  Solving  the  Unsteady  Three-Dimensional  Euler  and  Navier-Stokes 
Equations  on  Stationary  or  Dynamic  Grids",  Engineering  and  Industrial  Research  Station 
Repot*  ‘Mississippi  State  University,  No.  MSSU-EIRS-ASE-88-2,  February  1988. 

21 .  Arabshani,  A.  and  Whitfield,  D.  L.,  “A  Multiblock  Approach  to  Solving  the  Three-Dimensional 
Unsteady  Euler  Equations  about  A  Wing-Pylon-Store  Configuration:,  AIAA-89-3401-CP, 
August  1989. 

22.  Roe,  P.  L.,  “Approximate  Riemann  Solvers,  Parameter  Vector,  and  Difference  Schemes”, 
Journal  of  Computational  Physics,  Vol.  43.  p.  357, 1 981 . 

23.  Osher,  S.  and  Charavarthy,  S.,  “Very  High  Order  Accurate  TVD  Schemes”,  ICASE  Report 
No.  84-44,  September  1984. 

24.  Steger,  J.  L.  and  Warming,  R.  R,  “FJux  Vector  Splitting  of  the  Inviscid  Gasdynamic  Equations 
with  Applications  to  Finite  Difference  Methods",  Journal  of  Computational  Physics,  Vol.  40, 
p.  263, 1981. 

25.  Baldwin,  B.  S.  and  Lomax,  H.,  “Thir>-Layer  Approximation  and  Algebraic  Model  for  Sepa¬ 
rated  Turbulent  Flows",  AIAA-78-257,  January  1978. 

26.  Cooper,  G.  R.  and  Sirbaugh,  J.  R.,  “The  PARC  Distinction:  A  Practical  Flow  Simulator”, 
AIAA-90-2002,  AIAA/ASME/SAE/ASEE  26th  Joint  Propulsion  Conference,  Orlando,  FL, 
July  1990. 

27.  Pulliam,  J.  H.,  “Euler  and  Thin  Layer  Navier-Stokes  Codes:  ARC2D,  ARC3D”,  Notes  for 
Computational  Fluid  Dynamics  User’s  Workshop,  The  University  of  Tennessee  Space  Insti¬ 
tute,  Tullahoma,  TN,  UTSI  Publication  E02-4005-023-84,  p.  15.1,  March  12-16, 1984,. 

28.  Patel,  N.  R.,  Sturek,  W.  B.  and  Smith,  G.  A.,  “Parallel  Computation  of  Supersonic  Flows  Us¬ 
ing  a  Three-Dimensional,  Zonal,  Navier-Stokes  Code”,  BRL-TR-30XX,  November  1988. 

29.  Swanson,  R.  C.  andTurkel,  E.,  “Artificial  Dissipation  and  Central  Difference  Schemes  for  The 
Euler  and  Navier-Stokes  Equations”,  AIAA  8th  Computational  Fluid  Dynamics  Conference, 
Honolulu,  HA,  June  1987. 

30.  Soni,  B.  K.,  Weatherill,  N.  R,  Thompson,  J.  R,  “Grid  Adaptive  Strategies  in  CRD”,  Internation¬ 
al  Conference  on  HydroScience  &  Engineering,  Washington,  DC,  June  1 993. 

31 .  Soni,  B.  K.,  Huddleston,  D.  H.,  Arabshahi,  A.,  Vu,  B.,  Patel,  N.,  Clarke,  and  Coleman,  M.,  “A 
Study  of  CRD  Algorithms  Applied  to  Complete  Aircraft  Configurations”,  AIAA-93-0784, 
AIAA  31st  Aerospace  Sciences  Meeting,  Reno,  NV,  January  1993. 

32.  Yang,  J.  C.,  Soni,  B.  K.,  “Structured  Adaptive-Grid  Generation”,  accepted  for  publication. 
Journal  of  Applied  Mathematics  &  Computations,  1993. 

33.  Thompson,  J.  F.  and  Gatlin,  B.,  “Easy  EAGLE:  An  Introduction  to  the  EAGLE  Grid  Code”, 
Mississippi  State  University,  1988. 

34.  Thompson,  J.  R,  “A  General  Three-Dimensional  Elliptic  Grid  Generation  System  on  A  Com¬ 
posite  Block  Structure”,  Con. r  uter  Methods  in  Applied  Mechanics  &  Engineering,  p.  377, 
1987. 

35.  Remotigue,  M.  G.,  Hart,  E.  T.  and  Stokes,  M.  L.,  “EAGLEView:  A  Surface  and  Grid  Genera¬ 
tion  Program  and  Its  Data  Managemenf,  NASA,  CP-3143,  p.  319,  Proceedings  of  a  Work¬ 
shop,  NASA  Langley  Research  Center,  Hampton,  VA,  April  1992. 

36.  Soni,  B.,  “Grid  Generation  for  Internal  Flow  Configurations”,  Journal  of  Computers  &  Mathe¬ 
matics  with  Applications,  Vol.  25,  No.  5/6,  September  1992.  tt~ 

37.  Steinbrenner,  J.  P.  and  Chauner,  J.  R.,  “Recent  Enhancements  to  the  GRIDGEN  Structured 
Grid  Generation  System”,  NASA,  CP-3143,  p.  253,  Proceedings  of  a  Workshop,  NASA 
Langley  Research  Center,  Hampton,  VA,  April  1992. 

38.  Sorenson,  R.  L.,  “GRAPEVINE:  Grids  About  Anything  by  Poisson’s  Equations  in  A  Visually 
Interactive  Networking  Environment,  NASA,  CP-3143,  p.  319,  Proceedings  of  a  Workshop, 
NASA  Langley  Research  Center,  Hampton,  VA,  April  1992. 


29 


39.  Akdag,  V.  and  Wulf,  A.,  “Integrated  Geometry  and  Grid  Generation  System  for  Complex  Con¬ 
figurations”,  NASA  CP-3143,  p.  161 ,  Proceedings  of  a  Workshop,  Hampton,  VA,  April  1992. 

40.  Dener,  C.  and  Hirsch,  C.,  “IGG  -  An  Advanced  Interactive  Grid  Generation  System”,  p.  503, 
Proceedings  of  the  Third  International  Conference  on  Numerical  Grid  Generation  in  Com¬ 
putational  Fluid  Dynamics  and  Related  Fields.  Barcelona,  Spain,  June  1 991 . 

41 .  Soni,  B.  K.  and  Yang,  J.  L.,  “General  Purpose  Adaptive  Grid  System”,  AIAA-92-0664,  AIAA 
30th  Aerospace  Sciences  Meeting,  Reno,  NV,  January  1992. 

42.  Jameson,  A.,  Schmidt,  N.  and  Turkel,  E.,  “Numerical  Solutions  of  the  Euler  Equations  by  Fi¬ 
nite  Volume  Methods  Using  Runge-Kutta  Time-Stepping  Schemer",  AIAA-81-1259,  AIAA 
14th  Fluid  and  Plasma  Dynamics  Conference.  Palo  Alto.  CA,  1 981 . 

43.  Patel,  N.  R.,  Sturek,  W.  B.  and  Smith,  G.  A.,  “Parallel  Computation  of  Supersonic  Flows  Using 
a  Three-Dimensional,  Zonal,  Navier-Stokes  Code”,  BRL-TR-30XX,  November  1988. 

44.  Swanson,  R.  C.  and  Turkel,  E.,  “Artificial  Dissipation  and  Central  Diffemece  Schemes  for  the 
Euler  and  Navier-Stokes  Equations”,  AIAA  6th  Computational  Fluid  Dynamics  Conference, 
Honolulu,  HA,  June  1.987. 

45.  Soni,  B.  K.,  Weatherill,  N.  P.,  Thompson,  J.  F.,  “Structured-Unstructured  and  Adaptive  Grid 
Techniques  in  the  National  Grid  Projecf,  The  Second  U.  S.  Congress  on  Computational  Me¬ 
chanics,  June  1993. 

46.  Yu,  T.-Y,  “IGES  Transformer  and  NURBS  in  Grid  Generation”,  Master  Thesis  in  Aerosapce 
Engineering,  Mississippi  State  University,  Mississippi  State  University,  Mississippi,  MS 

1992. 

47.  Yoon,  Y.-H.,  “Enhancements  and  Extensions  of  EAGLE  Grid  Generation  System”,  Ph.  D. 
Dissertation,  Mississippi  State  University,  Mississippi  State,  MS,  1991. 

48.  Thompson,  J.  F.,  Warsi,  2.  U.  A.,  and  Mastin,  C.,  W.,  “Numerical  Grid  Generation:  Founda¬ 
tions  and  Applications”,  North  Holland  Publisher. 

49.  Soni,  B.  K.,  “Grid  Optimization:  A  Mixed  Approach”,  Proceedings  of  the  3rd  International 
Conference  of  Numerical  Grid  Generation  in  Computational  Fluid  Dynamics,  Barcelona, 
Spain,  June  1 991 ,  edited  by  A.  S.  Arcilla,  J.  Hauser,  P.  R.  Eiseman  and  J.  F.  Thompson:  North 
Holland. 

50.  Thompson,  J.  F.,  “A  Survey  of  Dynamically-Adaptive  Grids  in  the  Numerical  Solution  of  Par¬ 
tial  Differential  Equations”,  Applied  Numerical  Mathematics  1,  p.  3,  North-Holland,  1985. 

51 .  Connett,  W.  C.,  Aganwal,  R.  K.,  and  Schwartz,  A.  L.,  “An  Adaptive  Gnd-Generation  Scheme 
for  Flow  Field  Calculations”,  AIAA-87-0199. 

52.  Abolhassani,  J.  S.,  Smith,  R.  E.  and  Surendra,  N.  T.,  “Grid  Adaptation  for  Hypersonic  Flow”, 
AIAA-87-1169. 

53.  Yamaguchi,  R,  “Curves  and  Surfaces  in  Computer  Aided  Geometric  Design”,  Springer-Ver- 
lag,  1988. 

54.  Dwyer,  H.  A.,  “Grid  Adaptation  for  Problem  in  Fluid  Dynamics",  AIAA  Journal,  Vol.  22,  No. 
12,  December  1 984. 

55.  Eiseman,  P.  R.,  “Alternating  Direction  Adaptive  Grid  Generation”,  AIAA-83-1937, 1983. 

56.  Jameson,  A.,  Baker,  T.  J.  and  Weatherill,  N.  R,  “Calculation  of  Inviscid  Trat^nic  Flow  Over 
A  Complete  Aircraft”,  AIAA-86-01 03, 1 986. 

57.  Marchant,  M.  J.  and  Weatherill,  N.  R,  “The  Construction  of  Nearly  Orthogonal  Multiblock 
Grids  for  Compressible  Flow  Simulation”,  Communications  in  Applied  Numerical  Methods, 

1 993. 

58.  Nakashashi,  K.  and  Deiwert,  G.  S.,  ‘Three-Dimensional  Adaptive  Grid  Method”,  AIAA  Jour¬ 
nal,  Vol.  24,  p.  948,  1986. 


59.  Thompson,  J.  F.  and  Mastin,  C.  W,,  “Order  of  Difference  Expressions  in  Curvilinear  Coordi¬ 
nate  Systems",  Journal  of  Fluids  Engineering,  Vol.  107,  P.  241 , 1985. 

60.  Weatherill,  N.  P.  and  Soni,  B.,  “Grid  Generation  -  Refinement  in  Structured-Unstructured  Al¬ 
gorithms”,  Proceedings  of  the  3rd  International  Conference  of  Numerical  Grid  Generation  in 
Computational  Fluid  Dynamics,  Barcelona,  Spain,  June,  edited  by  A.  S.  Arcilla,  J.  Hauser, 
P.  R.  Eiseman  and  J.  R  Thompson,  North  Holland,  p.  143, 1991. 

61.  Weatherill,  N.  R,  “Mixed  Structured-Unstructured  Meshes  for  Aerodynamic  Flow  Simula¬ 
tion”,  The  Aeronautical  Journal,  Vol.  94,  p.  111-123,  No.  934, 1990. 

62.  Weatherill,  N.  R,  Hassan,  O.  and  Marcum,  D.  L.,  “Adaptive  Inviscid  Flow  Solutions  for  Aero¬ 
space  Geometries  on  Efficiently  Genreated  Unstructured  Tetrahedral  Meshes”, 
AIAA-93-0341 ,  AIAA  Aerospace  Sciences  Conference,  January  1993. 

63.  Kutler,  R,  “A  Perspective  of  Theoretical  and  Applied  Computational  Fluid  Dynamics”, 
AIAA-8S-0037,  1983. 


64.  Phares,  W.  J.,  Cooper,  G.  K.,  Swafford,  T.  W.,  Jones,  R.  R.,  “Application  of  Computational 
Fluid  Dynamics  to  Test  Facility  and  Experimental  Design”,  AIAA-86-1 733,  June  1986. 


65.  Droegemeier,  K.  K.,  “Numerical  Simulation  of  Thunderstorm  Outflows  and  Microbursts”, 
CRAY  CHANNELS,  p.  18,  Summer  1987. 

66.  Kent,  R.,  IGES  Editor,  IGES/PDES  Organization  Spline  Curves  &  Surfaces  in  the  Initial 
Graphics  Exchange  Specification  (IGES)  version  5.0  Appendix  B,  p.  417-422. 

67.  Soni,  B.  K.,  “Elliptic  Grid  Generation  System:  Control  Functions  Revisited”,  Journal  of  Ap¬ 
plied  Mathematics  &  Computers. 

68.  Soni,  B.  K.,  “An  Efficient  and  Robust  Algorithm  for  Grid  Refinemenr,  to  be  published. 

69.  Cooper,  G.  K.,  Sirbaugh,  J.  R.,  “PARC  Code:  Theory  and  Usage”,  AEDC-TR-89-1 5,  Arnold 
Engineering  Development  Center,  Arnold  AFB. 

70.  Yang,  J.  C.,  “General  Purpose  Adaptive  Grid  Generation  System”,  Ph.  D.  Dissertation,  Mis¬ 
sissippi  State  University,  Mississippi  State,  MS,  July  1993. 

71 .  Brackbill,  J.  U.  and  Saltzman,  J.  S.,  “Adaptive  Zoning  for  Singular  Problems  in  Two  Dimen¬ 
sion”,  Journal  of  Computational  Physics,  Vol.  46,  p.  342,  1 982. 


!•  ..»i;  I  !  •. 


.  j. Kirks; iiUfit. 

tire':  31 

my- 

'Ami 


vAv\\V5^vt',.j 


i",'. 


Figure  lb.  J  =  constant  grid  surface 


Figure  Ic.  I  =  constant  surface 
(I  =  173) 
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Figure  Id.  Grid  behavior  at  I  =  constant  cutting  through  wing 
Figure  1.  F15  Wing/Fuselage  Grid  (173  x  42  x  80) 


Figure  2a.  F15e  Vertical -Horizontal  Fin/Wing/Fuselage 


Figure  2b.  Cross  Sectional  View  of  Grid  and  Configuration 


Figure  3,  Vertical-Horizontal/Wing/Fuselage  Configurations 


ORIGINAL  GEOMETRY  DATA 


Figure  5.  Frontal  Profile  of  Original  CAD  Geometry 


Figure  6.  5-Block  Wing/Body/Tail  Configuration  With  Grid  Planes 


Figure  7.  Wing/Body/Tail  Configuration  Solution  As  Obtained  By  Ubiflow 
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GREEN- -OUTER INLET/N02ZLE  SURFACE 
PINK--LOWER  INLET/NOZZLE  SURFACF 
RED--UPPER  INLET/NOZZLE  SURFACF 
LIGHT  BLUE-  -INNET  INLET  SURFACE 


FRONTAL  VmWOF  SURFACE  GRID 

Figure  9.  Full  Surface  Grid  of  Body/Inlet/Nozzle  Configuration  Created  With  Genie-H- 


BLACK--FUSELAGE 

GREEN--OUTER  INLET/NOZZLE  SURFACE 
PINK-  -  LOWER  INLET/NOZZLE  SURFACE 
RED-  -  UPPER  INLET/NOZZLE  SURFACE 
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REAR  VIEW  OF  SURFACE  GRID 

Figure  10.  Rear  Perspective  of  Genie  Surface  Grid  Showing 
Nozzle  Geometry 
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Figure  12.  Detailed  Look  At  The  Inlet  Geometry  Modeled  Using  Genie++ 
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DETAIL  OF  NOZZLE  GRID 

Figure  13.  View  of  Nozzle  Surfaces 


Figure  15.  View  of  Overall  Volume  Grid  Showing  General  Point  Distribution  and  Outer  Boundary 
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Figure  18.  Grid  Plane  On  Interior  of  Inlet  to  Allow  for  Flow-Through  Calculations 


Figure  19.  I-Constant  Plane  Depicting  Grid  Lines  Entering  Inlet 
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Figure  24.  Arbitrary  Configuration  (with  close-up  grid 


Figure  25a.  Algebraic  Grid  Figure  25b.  Elliptic  Grid 
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Figure  27b.  Redistributed  Surface  Figure  27c.  Redistributed  Surface 

Cubic  Spline  Definition  Grid  NURBS  Definition 


Figure  28a.  C-Grid  on  A  Surface  Figure  29a.  Parametric  Space 

Associated  with  C-Typed  Grid 


Figure  28b.  Remapped  Q-Type  Grid 


Figure  29b. PBasame trie  Space  for 
0-Type  Grid 


Figure  30a.  Initial  0-Typ.e  Grid 


Figure  30c.  Remapped  Parametric 
Space  for  the  desired  H-Type  Grid 


Figure  30b.  Parametric  Space 
Associated  with  Initial  0-TypeviGrid 


Figure  30d.  Resulting  H-Type  Grid 
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Figure  39a.  Physical  Surface  Grid  Figure  39b.  Surface  Parametric  Space 
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Figure  41  .  Listing  of  Crossings  Found 
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.  Surface  Mach  Number  Comparison  Between  Roe  Algorithm  and 
Central  Difference  Algorithm 
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Figure  44.  Pressure  Comparison  on  Lower  Wing  Surface,  25%  from  Root 
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Figure  46.  Surface  Pressure  -  Parallel  Computation 


